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Claims 



1. A method for producing three-dimensional 
information of an object (4) in medical X-ray 
imaging, characterized in that 

- the object is X-radiated from at least two different 
directions and the said X-radiation is detected to form 
projection data of the object (4) 

- the object is modelled mathematically utilizing the 
projection data to solve the imaging geometry 
and/or the motion of the object, where the solving 
concerns either some or all parts of the imaging 
geometry and/or the motion of the object. 

- and said projection data and said mathematical 
modelling of the object are utilized in Bayesian 
inversion based on Bayes' formula 

Pirn) 

to produce three-dimensional information of the 
object the prior distribution Ppr(8) representing the 
prior knowledge of the imaging geometry and/or the 
motion of ihe object (4), the prior distribution p pr (x) 
representing mathematical modelling of the object, 
x representing the object image vector, which 
comprises values of the X-ray attenuation coefficient 
inside the object, 9 representing the parameter 
vector of the imaging geometry and/or the motion of . 
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the object (4), m representing projection data, the 
likelihood distribution p(mlx,e ) representing the X- 
radiation attenuation model between the object 
image vector x, geometry parameter vector 0, and 
projection data m, p(m) being a normalization 
constant and the posteriori distribution p(x,e|m) 
representing the three-dimensional information of 
the object (4) and the imaging geometry including 
the motion of the object. 

2. A method according to claim 1, characterized in that 
the three-dimensional information of the object (4) is 
one or more two-dimensional images representing x- 
ray attenuation coefficient along slices through the 
object. 

3- A method according to claim 1, characterized in that 
the three-dimensional information of the object (4) is a 
three-dimensional voxel representation of the X-ray 
attenuation in the object. 

4. A method according to claim 1, characterized in that 
the measurement model is m = A©x+e, where matrix Ae 
contains the weights describing how much each voxel 
contributes to the X-ray attenuation along the X-ray 
paths and the noise e is independent of object image 
vector x and the geometry parameter vector © leading 
to the likelihood distribution 

Pirn | x, 0) = p Misc (m - A^x) . 
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5. A method accordingjto claim 1, characterized in that 
the said mathematical modelling comprises that X- 
radiation attenuates when passing the object (4), which 
means that every image voxel is nonnegative. 

6. A method according! to claim 1, characterized in that 
mathematical modellinjg is expressed by the formula: 

Pp*(®) = exp(-c^y^C/iv(g?)) 

N 

where the sum is taken over a collection of 3D 
neighbourhoods N andl the value u N (*) depends only 

on the values of voxels belonging to the neighborhood 
N, and a is a positive r.egularization parameter used to 
tune the width of the prior distribution- 

i 

7- A method according to claim 1, characterized fn that 
the 3D tomographic problem is divided into a stack of 
2D tomographic probkjms and on every such 2D 
problem, the mathematical modelling is expressed by 
the formula: p ^ x) = eX p(-a ]T) U N (x)) 

N 



where the sum is taken over a collection of 2D 

i 

neighbourhoods N andl the value u N (x) depends only 

on the values of pixels belonging to the neighborhood 
N, and a is a positive regularization parameter used to 
tune the width of the prior distribution, and the 2D 
tomographic problems are related to each other by the 
formula 
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pr3D(x(j)) = exp(-YZSIxa)[k,q]-x(j-l)[k /q ] |), 
where the sums are taken over all pixels (k=l,...,K, 
q=l,... ,Q ) and y > 0 is another regularization 
parameter. 



8. A method according to claim 7, characterized in that 

i 

the neighborhoods consist of two adjacent pixels and U 
calculates a power of the absolute value of the 
difference, leading to the formula 

p pr {*™) = exp ( - «( jh |«<i>[A, g] - + 1, q ] |* - 

Aj=i q=l 
; K Q-l 

where s is a positive real number. 

9. A method according to claim 8, characterized in that 
s=l corresponding to total variation (TV) prior 
describing objects (4) consisting of different regions 
with well-defined boundaries. 

10. A method according to claim 1, characterized in 
that mathematical modelling is qualitative structural 
information of the target where the structural 
information is encoded, in prior distributions that are 
concentrated around object image vectors x that 
correspond to the physiological structures of the object 
(4)- 



11. A method according to claim 1, characterized in 
that mathematical modelling consists of a list or 
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probability distribution 1 of possible attenuation 
coefficient values in the object (4). 

12. A method according to claim 1/ characterized in 
that the X-ray imagingj geometry, such as X-ray source 
position, has unknown error modelled in the distribution 
p(m|x,9). 

13- A method according to claim 1, characterized in 
that the X-radiation measurement noise is Poisson 
distributed. 

i 

14. A method according to claim 1, characterized in 

i 

that the medical X-ray imaging is dental radiography. 

15- A method according to claim 1, characterized in 
that the medical X-ray imaging is surgical C-arm 
imaging. — — 

16- A method according to claim 1, characterized in 
that the medical x-ray imaging is mammography. 

s 

17. A method according to claim 1, characterized in 
that three-dimensional (information of the object (4) is 
produced on the basis of the maximum a posteriori 
estimator (MAP) which is calculated by the equation: 

P<3W I m > = max/>(;y I m ) 
m representing projection data and y=(x,9) 
representing the vector composed of the object image 
vector and the geometry parameter vector and where 
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i 

the maximum on the right hand side of the equation is 
taken over all y. 

18. A method according to claim 1, characterized in 

i 

that three-dimensional- information of the object (4) is 
produced on the basis 'of the conditional mean 
estimator (CM), which lis calculated by the equation: 

! 

where m represents projection data and y=(x,6) 
represents the vector composed of the object image 
vector and the geometry parameter vector. 

19. A medical X-ray device (5) arrangement for 
producing three-dimensional information of an object 
(4) in a medical X-ray Imaging, characterized in that the 
medical X-ray device (5) arrangement comprises: 

- an X-ray source (2) for X-radiating the object from 
at least two different directions 

- a detector (6) for detecting the X-radiation to form 
projection data of the object (4) 

- means (15) for modelling the object (4) 
mathematically utilizing the projection data 

to solve the imaging geometry and/or the motion of the 
object, where the solving concerns either some or all 
parts of the imaging geometry and/or the motion of the 
object. 

- and the medical X-ray device (5) arrangement 
comprises means (15) for utilizing said projection data 
and said mathematical -modelling of the object in 
Bayesian inversion based on Bayes' formula 
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p(ji g, m)= P,rfflP, f (*)P(H*.*) 

p{m) 

to produce three-dimensional information of the object, 
the prior distribution Pp r (e) representing the prior 
knowledge of the imaging geometry and/or the motion 
of the object (4), the prior distribution p Dr (x) 
representing mathematical modelling of the object, x 
representing the object image vector, which comprises 
values of the X-ray attenuation coefficient inside the 
object, 9 representing the parameter vector of the 
imaging geometry and/or the motion of the object (4), m 
representing projection data, the likelihood distribution 
p(m|x,e ) representing 1 the x-radiation attenuation 
model between the object image vector x, geometry 
parameter vector 0, and projection data m, p(m) being 
a normalization constant and the posteriori distribution 
p(x,6|m) representing the three-dimensional 
information of the object (4) and the imaging geometry 
including the motion of the object. 

20. A medical x-ray device (5) arrangement according 
to claim 19, characterized in that the three-dimensional 
information of the object (4) is one or more two- 
dimensional images representing X-ray attenuation 

coefficient along slices through the object. 

i 

21. A medical x-ray deyice (5) arrangement according 
to claim 19, characterized in that the three-dimensional 
information of the object (4) is a three-dimensional 
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voxel representation of the X-ray attenuation in the 
object. | 

22. A medical X-ray device (5) arrangement according 
to claim 19, characterized in that the medical X-ray 
device arrangement comprises means (15) for 
modelling the measurement as 

m = Aox+e, 

where matrix Ae contains the weights describing how 
much each voxel contributes to the X-ray attenuation 

along the X-ray paths and the noise e is independent of 

... ^ 

object image vector x and the geometry parameter 
vector 9 leading to the! likelihood distribution 

23. A medical X-ray device (5) arrangement according 
to claim 19 characterised in that the medical X-ray 
device arrangement comprises means (15) for 
modelling the object (4) mathematically so that X- 
radiation attenuates when passing the object (4), which 
means that every image voxel is nonnegative. 

i 

24. A medical X-ray device (5) arrangement according 
to claim 19, characterized in that the medical X-ray 
device arrangement comprises means (15) for 
modelling the object (4) mathematically by the formula: 

p VT {x) = exp(-ay^lOv-QK)) 

N 

where the sum is takeiji over a collection of 3D 
neighbourhoods N and! the value u s (x) depends only 
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on the values of voxels belonging to the neighborhood 
N, and a is a positive regularization parameter used to 
tune the width of the prior distribution- 

25, A medical x-ray device (5) arrangement according 
to claim 19, characterized in that the 3D tomographic 
problem is divided into a stack of 2D tomographic 
problems and on every such 2D problem, and the 
medical X-ray device arrangement comprises means 
(15) for modelling thejobject (4) mathematically by the 
formula: 

i>prO) = exp(-c* 

where the sum is taken over a collection of 2D 
neighbourhoods N and the value (J N (jc) depends only 

on the values of pixels; belonging to the neighborhood 
N r and a is a positive regularization parameter used to 
tune the width of the prior distribution, and the 2D 
tomographic problems are related to each other by the 
formula 

pr3D(xQ)) = exp(- Y I x(j)[k,q] - x(j-l)[k,q] [), 

where the sums are taiken over all pixels (k=l,...,K, 
q«l,...,Q ) and y > 0 is another regularization 
parameter. 

26. A medical X-ray device (5) arrangement according 
to claim 25, characterized in that the neighborhood 
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systems consist of twd neighboring pixels xj, xk or 
voxels xj, xk and u N (x) calculates a power of the 

i 

p pe (a^) = exp(- <*( k a) [k,q) ~ + 

fc=l q=l 

+ it 15 q ] - * w « +i ir 

absolute value of the difference, leading to the formula 
where s is a positive real number and a Is a 
regularization parameter used to tune the width of the 
prior distribution. 

i 

27. A medical X-ray device (5) arrangement according 
to claim 26, characterized in that the medical X-ray 
device arrangement comprises means (15) for 
modelling the object (4) mathematically by setting s=l 
corresponding to total! variation (JV) prior describing 
objects consisting of different regions with well-defined 
boundaries. 

28- A medical X-ray device (5) arrangement according 
to claim 19, characterized in that the medical X-ray 
device arrangement comprises means (15) for 
modelling the object (5) mathematically by assuming 
that mathematical modelling is qualitative structural 
information of the target where the structural 
information is encoded in prior distributions that are 
concentrated around image vectors x that correspond 
to the physiological structures of the target. 



i 

i 
i 
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29. A medical X-ray device (5) arrangement according 
to claim 19, characterized in that the medical X-ray 
device arrangement comprises means (15) for 
modelling the object (4) mathematically by assuming 
that mathematical modelling consists of a list of 

possible attenuation coefficient values in the object. 

i 
i 

30. A medical X-ray device (5) arrangement according 
to claim 19, characterized in that the medical x-ray 
device arrangement comprises means (15) for 
modelling the object (fl) mathematically by assuming 
that the X-ray imaging geometry, such as X-ray source 
position, has unknown 1 error modelled in the distribution 
p(m|x,e). 

31. A medical X-ray device (5) arrangement according 
to claim 19, characterized In that the medical X-ray 
device arrangement comprises means (15) for 
modelling the object (4) mathematically by assuming 
that X-radiation measurement noise is Poisson 
distributed. 

32. A medical X-ray device (5) arrangement according 
to claim 19, characterized in that the medical X-ray 
imaging is dental radiography. 

33. A medical X-ray device (5) arrangement according 
to claim 19, characterized in that the medical X-ray 
imaging is surgical C-arm imaging. 
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34, A medical X-ray device (5) arrangement according 
to claim 19, characterized in that the medical X-ray 
imaging is mammography. 



35. A medical X-ray device (5) arrangement according 
to claim 19, characterized in that the medical X-ray 
device arrangement comprises means (15) for 
producing three-dimensional information of the object 
(4) on the basis of the maximum a posteriori estimator 
(MAP), which is calculated by the equation: 

£>(3W I m) = max p(y \ rri) 



m representing 
representing the vector 
vector and the 
the maximum on the 
taken over all y. 



geometry 



projection data and y=(x,9) 

composed of the object image 
parameter vector, and where 
right hand side of the equation is 



36. A medical X-ray device (5) arrangement according 
to claim 19, characterized in that the medical X-ray 
device arrangement comprises means (15) for 
producing three-dimeipsional information of the object 
(4) on the basis of thct conditional mean estimator 
(CM), which is calculated by the equation 

y C M =fyp(y\™)<*y 



where m represents 
represents the vector 
vector and the geomdtry 



projection data and y=(x,9) 
composed of the object image 
parameter vector. 
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